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Abstract 

The generation of large-scale zonal flows by small-scale electrostatic 
drift waves in electron temperature gradient (ETG) driven turbulence model 
is considered. The generation mechanism is based on the modulational 
instability of a finite amplitude monochromatic drift wave. The threshold 
and growth rate of the instability as well as the optimal spatial scale of 
zonal flow are obtained. 



1 Introduction 

It is now an established fact that zonal flows (ZFs) (i.e., azimuthally symmet- 
ric modes that depend only on the radial coordinate) play a crucial role in 
regulating the nonlinear evolution of drift-wave instabilities in tokamaks, and 
consequently, the level of turbulent transport [1,2]. In particular, it is known 
that triggering of the L-H transition in tokamaks is related to the emergence 
of ZFs in the poloidal direction which suppresses the fluctuations and build up 
a barrier to the turbulent transport. It is widely thought that zonal flows are 
generated by the modulational instability of the turbulent spectrum of electro- 
static drift-wave perturbations [3,4]. Under this, two possible regimes can be 
indicated: (a) when the spectrum is broad, integration over all wavenumbers 
yields a resonant instability; (b) when the drift wave spectrum is narrow, one 
can consider the instability of a monochromatic drift wave and the instability is 
of the modulational type. In the latter case, standard schemes of the four-wave 
interaction can be applied [5-7]. 

Ion temperature gradient (ITG) modes and trapped electron modes are gen- 
erally regarded as the main candidates to explain the anomalous transport, but, 
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recently, there has been the growth of interest in electron temperature gradi- 
ent driven turbulence, produced by the electron temperature gradient driven 
modes [8-11]. This is related to the fact that the study of interaction of ETG 
modes with large-scale motions like ZFs or streamers (radially elongated vortex- 
type structures) is important for understanding the electron thermal transport 
within an internal transport barriers, when ITG turbulence is suppressed by 
E x B shear flow [12]. Note, that there is an important difference between 
the ITG and ETG models: in the ETG model we have the Boltzmann ion re- 
sponse for both waves and zonal flows, while the electron response to zonal flow 
perturbations is hydrodynamic in the ITG model [1, 13]. 

In the present paper we consider the excitation of ZFs by a finite amplitude 
monochromatic drift wave in the framework of ETG turbulence model. The 
corresponding ETG mode is assumed to be stable, i. e. we consider the region 
below the marginal stability boundary. We derive a set of coupled equations 
describing the nonlinear interaction of drift ETG modes and ZFs and show that 
ZFs can be readily excited by the modulational instability. 

This paper is organized as follows. In Sec. [5] a set of normalized fluid 
equations describing ETG drift modes - ZF interaction is introduced. In Sec. 
[31 the nonlinear dispersion relation is obtained and the modulational instability 
growth rate behavior is analyzed. Sec. [H contains summary and conclusions. 



2 Model equations 

Assuming a slab two-dimensional geometry, charge quasineutrality and the adi- 
abatic ion responce, we consider the following simplified model describing cur- 
vature driven ETG turbulence in the inviscid limit (see, e.g. [9]): 

^(l-A ± )^+|-(^ + P)-{^A^} = 0, (1) 

where tp and P are the normalized electrostatic potential and plasma pressure 
respectively, and {A, B} = d x Ad y B — d y Ad x B is the Jakobian. In equations 
(Q} and ([2]), following the notations of Ref. [9], we have rescaled the variables 
as follows 

r = — n — , 
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x', y' and t' being the original physical coordinates (with x' the poloidal and y' 
the radial coordinate), 

_ PsVt _ PsVt _ PsVt 

L n Lib Lip 

where L n , Lb and L p are the background gradient scales for the density, mag- 
netic field and pressure respectively, p s is the ion gyroradius calculated at the 
electron temperature T e , and r — Ti/T e . 

In the linear limit, equations |T]) and @ give the dispersion relation for ETG 
modes 

(3) 

where the plus sign describes the drift waves dispersion, while the minus sign 
corresponds to convective cells. The linear stability condition then reads as 

1 - 4r (k 2 + 1) > 0. 

If r > 1/4, the linear stability condition never holds, and ETG mode is always 
unstable. For r < 1/4, the stable ETG modes are confined inside the region 
k 2 < l/(4r) — 1. Below we restrict our analysis to the case of linearly stable 
ETG drift waves. 
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3 Modulational instability of ETG drift waves 
and zonal flow generation 

Assuming that the zonal flow varies on much larger timescale than ETG drift 
waves do, the standard decomposition into fast (~, ETG drift wave related) 
and slow (."?., zonal flow related) motions can be performed. The perturbations 
of electrostatic potential tp and plasma pressure P are presented as a sum of 
fast and slow parts 

ip = (p + (p, P = P + P. 

Following the standard averaging procedure (recall that we consider the region 
below the marginal stability boundary, i.e. the subcritical turbulence) one gets 
the following set of equations 

^(l-Ax)£={&A ± £}, (4) 
-| (1 - Aj_) <p + ^ (tp + Pj = {£, A ± <^} + A^} , (5) 

p ={P,?}, (6) 
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When obtaining the system (g)) - ([7]), we have assumed that the mean field is 
described by the same model equations as the fluctuations, but is driven by the 
Reynolds stress. This follows from the assumption of adiabatic ion response 
for the flow as well as fluctuations, which is valid for the ETG modes. We 
describe the interaction between drft ETG modes and ZF in terms of a four- 
wave coupling scheme, i.e., each fluctuation is taken to be coherent. Then, the 
ETG drift waves are considered as a superposition of the pump wave (k,uik) 
and two sidebands (k±,u)±), i.e., 



P = Pn + P + +P- 



where 



tpa = fo exp (ikr — iujktj + c.c, 
<p± = tp± exp (ik±r — uo±t\ + c.c, 
Pq = Pq exp ^ifcr — iw/fcij + c.c, 
P± = P± exp I ik±r — ioj±t J + c.c. 



(8) 

(9) 
(10) 

(11) 
(12) 



Zonal flow related electrostatic potential and pressure are taken in the form 



tp = ifq exp (iqr — iflt) + c.c, 
P = P q exp (iqr — iflt) + c.c. 



(13) 
(14) 



where q — (q, 0) is the wave vector of the zonal flow, and the resonant conditions 
k± = k ± q and u± = Uk ± hold, where Wk is the ETG drift mode frequency, 
given by equation §5§ with the " + " sign in front of the square root. 
Substitution of © - (HU) into the system flU) - © gives: 
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UJ+P+ + k +y rip + = iqk y [ip q P - foPq] , 


(19) 
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(20) 



Combining equations (fT5j) - (l20p . one can calculate the amplitudes of the up- 
shifted and down-shifted satellites 
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^ L q <p* (k 2 -q 2 ) + ^ faPg - 
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where we have introduced the notation 



(<^P * - <p* P q ) 



S± = lo± (fc| + 1) - k± + 



(24) 



Using the relation Pq = —(rk y /wk) ifo and eliminating Pq, from equations (fTS"|) , 
(fT5|) and (|21[) - (pM)) one can get the nonlinear dispersion relation in the form 



where 
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Generally, ZF dispersion which follows from the equation (j2l)|) is a fourth-order 
in f2, which, in principle, allow one to treat it analytically. We present below 
an analysis of nonlinear dispersion relation for the case when the ETG drift 
mode does not have the fc-component in the direction of inhomogeneity, i.e. for 
k x = 0. We also treated the case k x ^ and found that inclusion of nonzero k x 
does not introduce any qualitative change to the growth rate behavior except 
for the decrease of it's magnitude. Besides that, it turns out that in all studied 
cases the pump wave with k x = ensures the largest growth rate in the explored 
region of linear stability. 

For the case k x — 0, the nonlinear dispersion relation is reduced to the 
biquadratic equation 

where 



c A tt A + c 2 n 2 + c = 0, (26) 



CO = 4fc V (k 2 - q 2 ) (k 2 +q 2 + l) Cft 
+2k 2 q 4 [k 2 - 2k (l + 2k 2 ) u k + (l + 5fc 2 + 4k 4 - 2 (q 2 + q 4 )) uj 2 ] <p 2 

-q 4 (l + q 2 )co 4 , 

c 2 = 2k 2 q 2 (k 2 +q 2 + l) (k 2 + 2(q 4 + q 2 ) + l) $ 
+ (q 2 + 1) [k 2 - 4kcu k (k 2 + q 2 + 1) 
+2 (k 2 +q 2 + l) (2 (k 2 + 1) + q 2 ) LU 2 ] , 
C4 = -(g 2 + l) (k 2 +q 2 + l)\ 

which yields the modulational instability growth rate 



, , -c 2 ± \Jc\ - 4c c 4 , , 
7 = max Im | ±\j | . (27) 

The analysis of the growth rate l[27|) shows that the modulational instability 
can be excited (7 > 0) only above some threshold in the input power ipg. This 
threshold is a function of the scales A; and q of ETG drift wave and zonal flow. 
The minimum threshold corresponds to the case q — > and is given by 

^ - -7 ^ -v (M) 

2 - 2(1 + 2k 2 )v ph + (1 + k 2 )(l + 4k 2 )v 2 ph j 

where v p h = Wk/k is the phase velocity of ETG drift wave. As the zonal flow 
wavenumber q approaches k, the instability threshold goes to infinity, as it is 
illustrated in Figure [TJ 

Present analysis is restricted to the case when the pump ETG drift wave, as 
well as both sidebands are linearly stable, i.e. for each value of r parameter, all 
considered scales k and q are confined in the (k, q) - plane inside the circle 

k 2 + q 2 = -^-l. (29) 
4r 
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Figure 1: Schematic plot of the modulational instability threshold dependence 
on the scale of the excited zonal flow. ip\ h is given by 



Thus, the modulational instability parameter range under consideration is lim- 
ited. Instability domain in (fc, g)-planc is outlined in Figure O It is clear that 
circle shrinks to the origin when r approaches 1/4. 

In the Figure [3j modulational instability growth rates are plotted as a func- 
tions of k and q scales of ETG drift wave and zonal flow respectively for different 
values of the parameter r and pump wave amplitude ifiQ. Dashed lines in the 
contour planes confine the considered instability domain (schematically outlined 
in Figure [2]). Outside the circle given by the equation (f!?5)) . the growth rate was 
artificially forced to zero due to linear stability requirements. It is seen that as 
the parameter r grows, the region of instability contracts and the growth rate 
magnitude goes down. In all studied cases, the maximum growth rate always 
corresponded to the (k, q) - pair from the circle (|29p so that for the fixed k there 
is an optimal (i.e. corresponding to the maximum growth rate) scale q of the 
generated ZFs. With the increase of a pump wave amplitude, the growth rate 
as well as the instability domain increases, however in the region of physically 
appropriate values of ip^, no qualitative changes to the growth rate behavior 
were found. 



4 Summary and conclusions 

In the present paper, we have considered the nonlinear interaction between drift 
waves and ZFs in subcritical ETG turbulence. We have derived a set of equa- 
tions describing the dynamics of nonlinearly coupled ETG modes and ZFs. As- 
suming that the spectrum of fluctuations is sufficiently narrow, we analyzed the 
obtained set of equations in terms of a four-wave coupling scheme and obtained 
the nonlinear dispersion relation. We have shown that this dispersion relation 
predicts the modulational instability of a finite amplitude monochromatic drift 
wave. Thus, ETG drift fluctuations can be destabilized by the four-wave in- 
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Figure 2: Modulational instability domain in the (k, g)-plane is confined by the 
dashed lines. 



teraction mechanism with simultaneous generation of ZFs. We have found the 
threshold of the modulational instability and the dependence of the instabil- 
ity growth rate on spatial scales of ETG drift waves and excited ZFs. For the 
fixed wavenumber of the ETG mode k, the growth rate always has a maximum 
which is achieved for some intermediate value of ZF wavenumber q. When the 
amplitude of the ETG pump wave increases (as well as when the parameter r 
decreases), the region of modulational instability widens towards small fc-scales. 
The present results thus demonstrate that ZFs in subcritical ETG turbulence 
can be excited by the modulational instability. 
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